<!DOCTYPE html>
<html>
<head>
    <meta charset="UTF-8">
    <meta name="viewport" content="width=device-width, initial-scale=1, maximum-scale=1">
    <meta name="author" content="esy">
    
    <meta name="description" content="esy">
    
    
    
    
    
    
    <title>VF基础知识 | ESY</title>
    <link href="https://esyyes.github.io" rel="prefetch" />

    
<link rel="stylesheet" href="/css/bootstrap.min.css">
<link rel="stylesheet" href="/css/aos.css">
<link rel="stylesheet" href="/css/style.css">

    
<script src="/js/jquery.min.js"></script>

    
<script src="/js/bootstrap.min.js"></script>

    
<script src="/js/aos.js"></script>

    
<script src="/js/highslide/highslide-full.min.js"></script>

    
<link rel="stylesheet" href="/js/highslide/highslide.css">

    <style type="text/css">
        @media (max-width: 768px) {
            body {
                background-color: #f0f0f0;
                background: url('/imgs/xsbg.gif');
                background-attachment: fixed;
            }
        }
    </style>
    
    <!--<script type="text/javascript">
      if (document.images) {
        var avatar = new Image();
        avatar.src = '/imgs/avatar.jpg'
        var previews = 'preview1.jpg,preview2.jpg,preview3.jpg,preview4.jpg'.split(',')
        var previewsPreLoad = []
        for(var i = 0; i < length; i++) {
          previewsPreLoad.push(new Image())
          previewsPreLoad[previewsPreLoad.length - 1].src = '/imgs/preview' + previews[i]
        }
      }
    </script>-->
<meta name="generator" content="Hexo 5.2.0"></head>
<body>
    <!-- 背景轮播图功能 -->
    <section class="hidden-xs">
    <ul class="cb-slideshow">
        <li><span>天若</span></li>
        <li><span>有情</span></li>
        <li><span>天亦老</span></li>
        <li><span>我为</span></li>
        <li><span>长者</span></li>
        <li><span>续一秒</span></li>
    </ul>
</section>
    <!-- 欧尼酱功能, 谁用谁知道 -->
    
    <div class="gal-menu gal-dropdown">
    <div class="circle" id="gal">
        <div class="ring">
            <a href="https://esyyes.github.io" class="menuItem" style="left: 50%; top: 15%;">首页</a>
            
            <a class="menuItem" style="left: 80.3109%; top: 32.5%;">下一页</a>
            
            <a href="/archives" class="menuItem" style="left: 80.3109%; top: 67.5%;">归档</a>
            <a href="/about" class="menuItem" style="left: 50%; top: 85%;">关于</a>
            <a href="/message" class="menuItem" style="left: 19.6891%; top: 67.5%;">留言板</a>

            
            <a class="menuItem" style="left: 19.6891%; top: 32.5%;">上一页</a>
            
        </div>
        <audio id="audio" src="/imgs/oni.mp3"></audio>
    </div>
</div>
    
    <header class="navbar navbar-inverse" id="gal-header">
    <div class="container">
        <div class="navbar-header">
            <button type="button" class="navbar-toggle collapsed"
                    data-toggle="collapse" data-target=".bs-navbar-collapse"
                    aria-expanded="false">
                <span class="fa fa-lg fa-reorder"></span>
            </button>
            <a href="https://esyyes.github.io">
                
                <style>
                    #gal-header .navbar-brand {
                        height: 54px;
                        line-height: 24px;
                        font-size: 28px;
                        opacity: 1;
                        background-color: rgba(0,0,0,0);
                        text-shadow: 0 0 5px #fff,0 0 10px #fff,0 0 15px #fff,0 0 20px #228DFF,0 0 35px #228DFF,0 0 40px #228DFF,0 0 50px #228DFF,0 0 75px #228DFF;
                    }
                </style>
                <!-- 这里使用文字(navbar_text or config.title) -->
                <div class="navbar-brand">ESY</div>
                
            </a>
        </div>
        <div class="collapse navbar-collapse bs-navbar-collapse">
            <ul class="nav navbar-nav" id="menu-gal">
                
                
                <li class="">
                    <a href="/">
                        <i class="fa fa-home"></i>首页
                    </a>
                </li>
                
                
                
                <li class="">
                    <a href="/archives">
                        <i class="fa fa-archive"></i>归档
                    </a>
                </li>
                
                
                
                
                <li class="dropdown">
                    <!-- TODO 添加hover dropdown效果 -->
                    <a href="#" class="dropdown-toggle" data-toggle="dropdown"
                       aria-haspopup="true" aria-expanded="false" data-hover="dropdown">
                        <i class="fa fa-list"></i>分类
                    </a>
                    <ul class="dropdown-menu">
                        
                        
                        <li>
                            <a href="/categories/py-study/">py_study</a>
                        </li>
                        
                        <li>
                            <a href="/categories/nlp/">nlp</a>
                        </li>
                        
                        <li>
                            <a href="/categories/Graduation-work/">Graduation work</a>
                        </li>
                        
                        <li>
                            <a href="/categories/work/">work</a>
                        </li>
                        
                        <li>
                            <a href="/categories/hexo/">hexo</a>
                        </li>
                        
                        <li>
                            <a href="/categories/hexo%E5%AE%8C%E5%96%84/">-hexo完善</a>
                        </li>
                        
                        
                        <li>
                            <a href="/categories">...</a>
                        </li>
                        
                        
                    </ul>
                </li>
                
                
                
                
                
                <li class="dropdown">
                    <!-- TODO 添加hover dropdown效果 -->
                    <a href="#" class="dropdown-toggle" data-toggle="dropdown"
                       aria-haspopup="true" aria-expanded="false" data-hover="dropdown">
                        <i class="fa fa-tags"></i>标签
                    </a>
                    <ul class="dropdown-menu">
                        
                        
                        <li>
                            <a href="/tags/py-study/">py_study</a>
                        </li>
                        
                        <li>
                            <a href="/tags/nlp/">nlp</a>
                        </li>
                        
                        <li>
                            <a href="/tags/Graduation-work/">Graduation work</a>
                        </li>
                        
                        <li>
                            <a href="/tags/work/">work</a>
                        </li>
                        
                        <li>
                            <a href="/tags/hexo/">hexo</a>
                        </li>
                        
                        <li>
                            <a href="/tags/%E4%B8%AA%E4%BA%BA%E5%8D%9A%E5%AE%A2%E6%90%AD%E5%BB%BA/">-个人博客搭建</a>
                        </li>
                        
                        
                        <li>
                            <a href="/tags">...</a>
                        </li>
                        
                        
                    </ul>
                </li>
                
                
                
                
                <li class="">
                    <a href="/about">
                        <i class="fa fa-user"></i>关于我
                    </a>
                </li>
                
                
            </ul>
        </div>
    </div>
</header>
    <div id="gal-body">
        <div class="container">
            <div class="row">
                <div class="col-md-8 gal-right" id="mainstay">
                    
<article class="article well article-body" id="article">
    <div class="breadcrumb">
        <i class="fa fa-home"></i>
        <a href="https://esyyes.github.io">ESY</a>
        >
        <span>VF基础知识</span>
    </div>
    <!-- 大型设备详细文章 -->
    <div class="hidden-xs">
        <div class="title-article">
            <h1>
                <a href="/2020/05/18/python%20work/VF%E5%9F%BA%E7%A1%80%E7%9F%A5%E8%AF%86/">VF基础知识</a>
            </h1>
        </div>
        <div class="tag-article">
            
            <span class="label label-gal">
                <i class="fa fa-calendar"></i> 2020-05-18
            </span>
            
        </div>
    </div>
    <!-- 小型设备详细文章 -->
    <div class="visible-xs">
        <center>
            <div class="title-article">
                <h4>
                    <a href="/2020/05/18/python%20work/VF%E5%9F%BA%E7%A1%80%E7%9F%A5%E8%AF%86/">VF基础知识</a>
                </h4>
            </div>
            <p>
                <i class="fa fa-calendar"></i> 2020-05-18
            </p>
            <p>
                
                
            </p>
        </center>
    </div>
    <div class="content-article">
        <h2 id="心室纤颤-Ventricular-Fibrillation-，VF"><a href="#心室纤颤-Ventricular-Fibrillation-，VF" class="headerlink" title="心室纤颤(Ventricular Fibrillation ，VF)"></a>心室纤颤(Ventricular Fibrillation ，VF)</h2><p>心室纤颤(Ventricular Fibrillation ，VF)特点：</p>
<p>QRS-T波群完全消失，代之以形态不同、大小各异间距极不匀齐的颤动波(f波)，频率为250～500次/分钟，颤动波之间无等电线。</p>
<p><strong>心室扑动典型的心电图特点</strong></p>
<p>连续而规则、宽大、畸形的QRS波，即心室扑动波。QRS波的时限长，在0.12s以上，QRS波呈向上向下的波幅似正弦样曲线与T波无法分开，QRS波之间无等电线。QRS波频率多在180～250次/分钟，有时可低达150次/分钟或高达300次/分钟。P波消失。</p>
<h2 id="心室纤颤的数据库："><a href="#心室纤颤的数据库：" class="headerlink" title="心室纤颤的数据库："></a>心室纤颤的数据库：</h2><p><img src="H:\myboke\mybike\source\images\心电数据库.png" alt="image-20200519141433962"></p>
<p><img src="H:\myboke\mybike\source\images\vfdb介绍.png" alt="image-20200519141519264"></p>
<h2 id="数据长短："><a href="#数据长短：" class="headerlink" title="数据长短："></a>数据长短：</h2><p><img src="H:\myboke\mybike\source\images\分析数据长短.png" alt="image-20200519141626942"></p>
<p><img src="H:\myboke\mybike\source\images\MITDB数据选择" alt="image-20200519141744095"></p>
<h2 id="总结："><a href="#总结：" class="headerlink" title="总结："></a>总结：</h2><p>直接选择数据库为VFDB，室颤的，然后根据标签选择室颤的信号，筛选10个的样子，每个10s，然后生成excel表格，再选择正常的数据10段，最后来区分，利用利用python直接来读取数据，然后生成数据组合!再进行预处理。最后筛选出一个算法进行进行分类，并有个准确率输出。</p>
<h1 id="VF数据库"><a href="#VF数据库" class="headerlink" title="VF数据库"></a>VF数据库</h1><p>含有室颤的数据有：424,426,430</p>
<p>数据长短为：2100s，35min</p>
<p>采样率为：250</p>
<p>首先确定VF和VFIB的位置，将其生成excel文档，10s长短，即2500数据点，</p>
<h2 id="python-wfdb截取数据"><a href="#python-wfdb截取数据" class="headerlink" title="python-wfdb截取数据"></a>python-wfdb截取数据</h2><figure class="highlight python"><table><tr><td class="gutter"><pre><span class="line">1</span><br><span class="line">2</span><br><span class="line">3</span><br><span class="line">4</span><br><span class="line">5</span><br><span class="line">6</span><br><span class="line">7</span><br><span class="line">8</span><br><span class="line">9</span><br><span class="line">10</span><br><span class="line">11</span><br><span class="line">12</span><br><span class="line">13</span><br><span class="line">14</span><br><span class="line">15</span><br><span class="line">16</span><br><span class="line">17</span><br><span class="line">18</span><br><span class="line">19</span><br><span class="line">20</span><br><span class="line">21</span><br><span class="line">22</span><br><span class="line">23</span><br><span class="line">24</span><br><span class="line">25</span><br><span class="line">26</span><br><span class="line">27</span><br><span class="line">28</span><br><span class="line">29</span><br><span class="line">30</span><br><span class="line">31</span><br><span class="line">32</span><br><span class="line">33</span><br><span class="line">34</span><br></pre></td><td class="code"><pre><span class="line"><span class="comment"># -*- coding: utf-8 -*-</span></span><br><span class="line"><span class="comment"># @Time     : 2020/5/19</span></span><br><span class="line"><span class="comment"># @Author   : esy</span></span><br><span class="line"></span><br><span class="line"><span class="keyword">import</span> wfdb <span class="keyword">as</span> wfdb</span><br><span class="line"><span class="keyword">import</span> pandas <span class="keyword">as</span> pd</span><br><span class="line"><span class="keyword">import</span> matplotlib.pyplot <span class="keyword">as</span> plt</span><br><span class="line"><span class="keyword">import</span> warnings</span><br><span class="line"></span><br><span class="line"><span class="comment"># 忽略警告</span></span><br><span class="line">warnings.filterwarnings(<span class="string">&quot;ignore&quot;</span>)</span><br><span class="line"></span><br><span class="line"><span class="comment"># 根据标签筛选出室颤信号</span></span><br><span class="line">record = wfdb.rdrecord(<span class="string">&#x27;F:/VF/VFDB/424&#x27;</span>, channels=[<span class="number">0</span>], sampfrom=<span class="number">314749</span>, sampto=<span class="number">341173</span>)</span><br><span class="line">VF = record.p_signal</span><br><span class="line"></span><br><span class="line">vf = []</span><br><span class="line"><span class="keyword">for</span> i <span class="keyword">in</span> <span class="built_in">range</span>(record.sig_len):</span><br><span class="line">    vf.append(VF[i][<span class="number">0</span>])</span><br><span class="line">vf = vf[<span class="number">5000</span>:<span class="number">7500</span>]</span><br><span class="line"></span><br><span class="line">fig = plt.figure(figsize=(<span class="number">7</span>, <span class="number">4</span>))</span><br><span class="line">plt.plot(vf, linewidth=<span class="number">1</span>)</span><br><span class="line"></span><br><span class="line"></span><br><span class="line">plt.ylabel(<span class="string">&quot;U/mv&quot;</span>, fontsize=<span class="number">10</span>)</span><br><span class="line">plt.xticks(fontsize=<span class="number">8</span>)</span><br><span class="line">plt.xlabel(<span class="string">&quot;time/s&quot;</span>, fontsize=<span class="number">10</span>)</span><br><span class="line">plt.yticks(fontsize=<span class="number">8</span>)</span><br><span class="line">plt.show()</span><br><span class="line"></span><br><span class="line"></span><br><span class="line">data = pd.DataFrame(vf)</span><br><span class="line">data.to_excel(<span class="string">&quot;vf1.xlsx&quot;</span>)</span><br></pre></td></tr></table></figure>

<p><img src="H:\myboke\mybike\source\images\VF的数据样本10s的】.png" alt="image-20200520104746150"></p>
<p><strong>筛选10个数据片段</strong></p>
<p><strong>筛选正常数据</strong></p>
<p>下一步就是用算法将其判断出来，写个循环看看，Ndata数据的效果，用pt算法，另外一个也试试</p>
<p>数据片段的筛选，重新保存筛选，先选择算法进行分析。</p>
<p>TCI算法：</p>
<p><a target="_blank" rel="noopener" href="https://link.springer.com/article/10.1186/1475-925X-4-60">https://link.springer.com/article/10.1186/1475-925X-4-60</a></p>
<p>随便设定一个大于250和小于250就行了</p>
<p> 纪念下人生第一单哈哈哈</p>
<figure class="highlight plain"><table><tr><td class="gutter"><pre><span class="line">1</span><br><span class="line">2</span><br><span class="line">3</span><br><span class="line">4</span><br><span class="line">5</span><br><span class="line">6</span><br><span class="line">7</span><br><span class="line">8</span><br><span class="line">9</span><br><span class="line">10</span><br><span class="line">11</span><br><span class="line">12</span><br><span class="line">13</span><br><span class="line">14</span><br><span class="line">15</span><br><span class="line">16</span><br><span class="line">17</span><br><span class="line">18</span><br><span class="line">19</span><br><span class="line">20</span><br><span class="line">21</span><br><span class="line">22</span><br><span class="line">23</span><br><span class="line">24</span><br><span class="line">25</span><br><span class="line">26</span><br><span class="line">27</span><br><span class="line">28</span><br><span class="line">29</span><br><span class="line">30</span><br><span class="line">31</span><br><span class="line">32</span><br><span class="line">33</span><br><span class="line">34</span><br><span class="line">35</span><br><span class="line">36</span><br><span class="line">37</span><br><span class="line">38</span><br><span class="line">39</span><br><span class="line">40</span><br><span class="line">41</span><br><span class="line">42</span><br><span class="line">43</span><br><span class="line">44</span><br></pre></td><td class="code"><pre><span class="line">clc;clear;close all</span><br><span class="line">% 读取室颤数据</span><br><span class="line">data &#x3D; xlsread(&#39;F:\VF\VFdata\0.xlsx&#39;);</span><br><span class="line">% 读取正常数据</span><br><span class="line">% data &#x3D; xlsread(&#39;F:\VF\Ndata\0.xlsx&#39;);</span><br><span class="line"></span><br><span class="line">% 采样频率</span><br><span class="line">fs &#x3D; 250;</span><br><span class="line"></span><br><span class="line">% 想要采样的时间长度为10s，所以数据长短通通为2500</span><br><span class="line">vf &#x3D; data(2:2501, 2);</span><br><span class="line"></span><br><span class="line">% 将采样点的横坐标转换为s</span><br><span class="line">N &#x3D; length(vf);</span><br><span class="line">m &#x3D; 0 : N-1;</span><br><span class="line">t &#x3D; m &#x2F; fs;</span><br><span class="line"></span><br><span class="line">% 绘制出原始数据</span><br><span class="line">subplot(2,1,1)</span><br><span class="line">plot(t,vf,&#39;b&#39;);</span><br><span class="line"></span><br><span class="line">title(&#39;原始信号&#39;);</span><br><span class="line">xlabel(&#39;time&#x2F;s&#39;);</span><br><span class="line">ylabel(&#39;U&#x2F;mV&#39;);</span><br><span class="line"></span><br><span class="line">% 可要可不要这一步将数据平滑一下</span><br><span class="line">subplot(2,1,2)</span><br><span class="line">df &#x3D; smooth(vf,3);</span><br><span class="line">plot(t,df);</span><br><span class="line">title(&#39;平滑后的信号&#39;);</span><br><span class="line">xlabel(&#39;time&#x2F;s&#39;);</span><br><span class="line">ylabel(&#39;U&#x2F;mV&#39;);</span><br><span class="line"></span><br><span class="line">% 利用TCI (Threshold crossing intervals)算法对TCI值进行判断</span><br><span class="line">% verbose &#x3D; 0数据直接输出8个3s窗口的平均值，verbose &#x3D; 1，进入调试界面</span><br><span class="line">verbose &#x3D; 1;</span><br><span class="line">tci &#x3D; VF_TCI(vf,fs,10,verbose);</span><br><span class="line"></span><br><span class="line">% VF的频率为250～500次&#x2F;分钟，转换为ms为120-240ms，即可判定为室颤</span><br><span class="line">if (tci &gt;&#x3D; 120) &amp;&amp; (tci &lt;&#x3D; 240)</span><br><span class="line">    disp(&#39;该信号为室颤信号&#39;)</span><br><span class="line">else</span><br><span class="line">    disp(&#39;该信号为正常信号&#39;)</span><br><span class="line">end</span><br></pre></td></tr></table></figure>

<p>读取单个的图</p>
<figure class="highlight matlab"><table><tr><td class="gutter"><pre><span class="line">1</span><br><span class="line">2</span><br><span class="line">3</span><br><span class="line">4</span><br><span class="line">5</span><br><span class="line">6</span><br><span class="line">7</span><br><span class="line">8</span><br><span class="line">9</span><br><span class="line">10</span><br><span class="line">11</span><br><span class="line">12</span><br><span class="line">13</span><br><span class="line">14</span><br><span class="line">15</span><br><span class="line">16</span><br><span class="line">17</span><br><span class="line">18</span><br><span class="line">19</span><br><span class="line">20</span><br><span class="line">21</span><br><span class="line">22</span><br><span class="line">23</span><br><span class="line">24</span><br><span class="line">25</span><br><span class="line">26</span><br><span class="line">27</span><br><span class="line">28</span><br><span class="line">29</span><br><span class="line">30</span><br><span class="line">31</span><br><span class="line">32</span><br><span class="line">33</span><br><span class="line">34</span><br><span class="line">35</span><br><span class="line">36</span><br><span class="line">37</span><br><span class="line">38</span><br><span class="line">39</span><br><span class="line">40</span><br><span class="line">41</span><br><span class="line">42</span><br><span class="line">43</span><br><span class="line">44</span><br><span class="line">45</span><br><span class="line">46</span><br><span class="line">47</span><br><span class="line">48</span><br></pre></td><td class="code"><pre><span class="line">clc;clear;close all</span><br><span class="line"></span><br><span class="line"><span class="comment">% 数据标签，将VF数据判定为1，正常的判定为0</span></span><br><span class="line">tag = [];</span><br><span class="line"><span class="keyword">for</span> <span class="built_in">i</span> = <span class="number">0</span> : <span class="number">31</span></span><br><span class="line">    <span class="keyword">if</span> <span class="built_in">i</span> &lt;= <span class="number">18</span></span><br><span class="line">        tag(<span class="built_in">i</span>+<span class="number">1</span>,<span class="number">1</span>) = <span class="number">0</span>;</span><br><span class="line">    <span class="keyword">else</span></span><br><span class="line">        tag(<span class="built_in">i</span>+<span class="number">1</span>,<span class="number">1</span>) = <span class="number">1</span>;</span><br><span class="line">    <span class="keyword">end</span></span><br><span class="line"><span class="keyword">end</span></span><br><span class="line"></span><br><span class="line"><span class="comment">% 循环读取数据，将VF数据判定为1，正常的判定为0</span></span><br><span class="line">num_vf = <span class="number">0</span>;</span><br><span class="line">num_n = <span class="number">0</span>;</span><br><span class="line">all = [];</span><br><span class="line">TCI = [];</span><br><span class="line"><span class="keyword">for</span> <span class="built_in">i</span> = <span class="number">0</span> : <span class="number">31</span></span><br><span class="line">    data = xlsread([<span class="string">&#x27;F:\VF\all_data\&#x27;, num2str(i), &#x27;</span>.xlsx&#x27;]);</span><br><span class="line">    fs = <span class="number">250</span>;</span><br><span class="line">    vf = data(<span class="number">2</span>:<span class="number">2501</span>, <span class="number">2</span>);    </span><br><span class="line">    tci = VF_TCI(vf,fs,<span class="number">10</span>,<span class="number">0</span>);</span><br><span class="line">    TCI(<span class="built_in">i</span>+<span class="number">1</span>,<span class="number">1</span>) = tci;</span><br><span class="line">    <span class="keyword">if</span> (tci &gt;= <span class="number">120</span>) &amp;&amp; (tci &lt;= <span class="number">240</span>)</span><br><span class="line">        num_vf = num_vf + <span class="number">1</span>;</span><br><span class="line">        all(<span class="built_in">i</span>+<span class="number">1</span>,<span class="number">1</span>) = <span class="number">1</span>;</span><br><span class="line">    <span class="keyword">else</span></span><br><span class="line">        num_n = num_n + <span class="number">1</span>;</span><br><span class="line">        all(<span class="built_in">i</span>+<span class="number">1</span>,<span class="number">1</span>) = <span class="number">0</span>;</span><br><span class="line">    <span class="keyword">end</span></span><br><span class="line"><span class="keyword">end</span></span><br><span class="line"></span><br><span class="line"><span class="comment">% 进行比较，用标签和这个你判定的数据进行比较</span></span><br><span class="line">ALL = <span class="number">32</span>;</span><br><span class="line">True = <span class="number">0</span>;</span><br><span class="line">F = <span class="number">0</span>;</span><br><span class="line"><span class="keyword">for</span> <span class="built_in">i</span> = <span class="number">0</span> : <span class="number">31</span></span><br><span class="line">    <span class="keyword">if</span> all(<span class="built_in">i</span>+<span class="number">1</span>,<span class="number">1</span>) == tag(<span class="built_in">i</span>+<span class="number">1</span>,<span class="number">1</span>)</span><br><span class="line">        True = True + <span class="number">1</span>;</span><br><span class="line">    <span class="keyword">else</span></span><br><span class="line">        F = F + <span class="number">1</span>;</span><br><span class="line">    <span class="keyword">end</span>     </span><br><span class="line"><span class="keyword">end</span></span><br><span class="line"></span><br><span class="line"><span class="comment">% 输出得到判定准确的值，差不多准确率就基本可以了，毕竟100%了，再转换为百分数</span></span><br><span class="line">ACR = True / ALL;</span><br><span class="line">accuracy = strcat(num2str(ACR * <span class="number">100</span>),<span class="string">&#x27;%&#x27;</span>);</span><br><span class="line"><span class="built_in">disp</span>([<span class="string">&#x27;准确率为:&#x27;</span>,accuracy]);</span><br></pre></td></tr></table></figure>

<p>输出为准确率</p>
<figure class="highlight matlab"><table><tr><td class="gutter"><pre><span class="line">1</span><br><span class="line">2</span><br><span class="line">3</span><br><span class="line">4</span><br><span class="line">5</span><br><span class="line">6</span><br><span class="line">7</span><br><span class="line">8</span><br><span class="line">9</span><br><span class="line">10</span><br><span class="line">11</span><br><span class="line">12</span><br><span class="line">13</span><br><span class="line">14</span><br><span class="line">15</span><br><span class="line">16</span><br><span class="line">17</span><br><span class="line">18</span><br><span class="line">19</span><br><span class="line">20</span><br><span class="line">21</span><br><span class="line">22</span><br><span class="line">23</span><br><span class="line">24</span><br><span class="line">25</span><br><span class="line">26</span><br><span class="line">27</span><br><span class="line">28</span><br><span class="line">29</span><br><span class="line">30</span><br><span class="line">31</span><br><span class="line">32</span><br><span class="line">33</span><br><span class="line">34</span><br><span class="line">35</span><br><span class="line">36</span><br><span class="line">37</span><br><span class="line">38</span><br><span class="line">39</span><br><span class="line">40</span><br><span class="line">41</span><br><span class="line">42</span><br><span class="line">43</span><br><span class="line">44</span><br><span class="line">45</span><br><span class="line">46</span><br><span class="line">47</span><br><span class="line">48</span><br><span class="line">49</span><br><span class="line">50</span><br><span class="line">51</span><br><span class="line">52</span><br><span class="line">53</span><br><span class="line">54</span><br><span class="line">55</span><br><span class="line">56</span><br><span class="line">57</span><br><span class="line">58</span><br><span class="line">59</span><br><span class="line">60</span><br><span class="line">61</span><br><span class="line">62</span><br><span class="line">63</span><br><span class="line">64</span><br><span class="line">65</span><br><span class="line">66</span><br><span class="line">67</span><br><span class="line">68</span><br><span class="line">69</span><br><span class="line">70</span><br><span class="line">71</span><br><span class="line">72</span><br><span class="line">73</span><br><span class="line">74</span><br><span class="line">75</span><br><span class="line">76</span><br><span class="line">77</span><br><span class="line">78</span><br><span class="line">79</span><br><span class="line">80</span><br><span class="line">81</span><br><span class="line">82</span><br><span class="line">83</span><br><span class="line">84</span><br><span class="line">85</span><br><span class="line">86</span><br><span class="line">87</span><br><span class="line">88</span><br><span class="line">89</span><br><span class="line">90</span><br><span class="line">91</span><br><span class="line">92</span><br><span class="line">93</span><br><span class="line">94</span><br><span class="line">95</span><br><span class="line">96</span><br><span class="line">97</span><br><span class="line">98</span><br><span class="line">99</span><br><span class="line">100</span><br><span class="line">101</span><br><span class="line">102</span><br><span class="line">103</span><br><span class="line">104</span><br><span class="line">105</span><br><span class="line">106</span><br><span class="line">107</span><br><span class="line">108</span><br><span class="line">109</span><br><span class="line">110</span><br><span class="line">111</span><br><span class="line">112</span><br><span class="line">113</span><br><span class="line">114</span><br><span class="line">115</span><br><span class="line">116</span><br><span class="line">117</span><br><span class="line">118</span><br><span class="line">119</span><br><span class="line">120</span><br><span class="line">121</span><br><span class="line">122</span><br><span class="line">123</span><br></pre></td><td class="code"><pre><span class="line"><span class="function"><span class="keyword">function</span> <span class="title">tci</span> = <span class="title">VF_TCI</span><span class="params">(xf,fs,wL,verbose)</span></span></span><br><span class="line"></span><br><span class="line"><span class="comment">% INPUT:</span></span><br><span class="line"><span class="comment">% - xf: ecg signal (preprocessed)</span></span><br><span class="line"><span class="comment">% - fs: sampling frequency</span></span><br><span class="line"><span class="comment">% - wL: window length, in seconds </span></span><br><span class="line"><span class="comment">% - verbose: debugging variable (1: plot; 0: default, not ploting)</span></span><br><span class="line"><span class="comment">%</span></span><br><span class="line"><span class="comment">% OUTPUT</span></span><br><span class="line"><span class="comment">% - tci parameter</span></span><br><span class="line"></span><br><span class="line">wl = <span class="number">1</span>*fs;  <span class="comment">% 1-sec window samples</span></span><br><span class="line">wa = <span class="number">3</span>*fs;  <span class="comment">% 3-sec window samples</span></span><br><span class="line"></span><br><span class="line"></span><br><span class="line">L = wL<span class="number">-3</span>+<span class="number">1</span>; <span class="comment">% number of 3-sec windows in wl segment </span></span><br><span class="line"></span><br><span class="line"><span class="comment">%becg = zeros(1,wa);</span></span><br><span class="line">becg1 = <span class="built_in">zeros</span>(<span class="number">1</span>,wl);</span><br><span class="line">becg2 = <span class="built_in">zeros</span>(<span class="number">1</span>,wl);</span><br><span class="line">becg3 = <span class="built_in">zeros</span>(<span class="number">1</span>,wl);</span><br><span class="line"></span><br><span class="line">tci6 = <span class="built_in">zeros</span>(<span class="number">1</span>,L);</span><br><span class="line"></span><br><span class="line"><span class="keyword">for</span> <span class="built_in">j</span>=<span class="number">0</span>:L<span class="number">-1</span></span><br><span class="line">    wsamples1 = <span class="built_in">j</span>*wl+<span class="number">1</span>:(<span class="built_in">j</span>+<span class="number">1</span>)*wl;</span><br><span class="line">    wsamples2 = (<span class="built_in">j</span>+<span class="number">1</span>)*wl+<span class="number">1</span>:(<span class="built_in">j</span>+<span class="number">2</span>)*wl;</span><br><span class="line">    wsamples3 = (<span class="built_in">j</span>+<span class="number">2</span>)*wl+<span class="number">1</span>:(<span class="built_in">j</span>+<span class="number">3</span>)*wl;</span><br><span class="line">    </span><br><span class="line">    stage1 = xf(wsamples1)-<span class="built_in">mean</span>(xf(wsamples1)); maxv = <span class="built_in">max</span>(stage1); </span><br><span class="line">    th1 = <span class="number">0.2</span>*maxv; becg1(stage1&gt;th1) = <span class="number">1</span>;</span><br><span class="line">    stage2 = xf(wsamples2)-<span class="built_in">mean</span>(xf(wsamples2)); maxv = <span class="built_in">max</span>(stage2); </span><br><span class="line">    th2 = <span class="number">0.2</span>*maxv; becg2(stage2&gt;th2) = <span class="number">1</span>;</span><br><span class="line">    stage3 = xf(wsamples3)-<span class="built_in">mean</span>(xf(wsamples3)); maxv = <span class="built_in">max</span>(stage3); </span><br><span class="line">    th3 = <span class="number">0.2</span>*maxv; becg3(stage3&gt;th3) = <span class="number">1</span>;</span><br><span class="line">    </span><br><span class="line">    becg = [becg1 becg2 becg3];</span><br><span class="line">    </span><br><span class="line">    aux = [<span class="number">0</span> diff(becg)];</span><br><span class="line">    </span><br><span class="line">    s1 = <span class="built_in">find</span>(aux(<span class="number">1</span>:wl)==<span class="number">-1</span>);</span><br><span class="line">    <span class="keyword">if</span> <span class="built_in">isempty</span>(s1)</span><br><span class="line">        t1 = <span class="number">1</span>;</span><br><span class="line">    <span class="keyword">else</span></span><br><span class="line">        t1 = (wl-s1(<span class="keyword">end</span>))/fs;</span><br><span class="line">    <span class="keyword">end</span></span><br><span class="line">    </span><br><span class="line">    index = <span class="built_in">find</span>(aux(wl+<span class="number">1</span>:<span class="number">2</span>*wl));</span><br><span class="line">    s2 = aux(wl+<span class="number">1</span>:<span class="number">2</span>*wl);</span><br><span class="line">    pulses = s2(index);</span><br><span class="line">    </span><br><span class="line">    <span class="keyword">if</span> pulses(<span class="number">1</span>) == <span class="number">-1</span> &amp;&amp; pulses(<span class="keyword">end</span>) == <span class="number">-1</span></span><br><span class="line">        t2 = <span class="number">0</span>;</span><br><span class="line">        t3 = (wl-index(<span class="keyword">end</span>))/fs;</span><br><span class="line">        N = (<span class="built_in">length</span>(pulses)+<span class="number">1</span>)/<span class="number">2</span>;</span><br><span class="line"></span><br><span class="line">    <span class="keyword">elseif</span> pulses(<span class="number">1</span>) == <span class="number">1</span> &amp;&amp; pulses(<span class="keyword">end</span>) == <span class="number">1</span></span><br><span class="line">        t2 = index(<span class="number">1</span>)/fs;</span><br><span class="line">        t3 = <span class="number">0</span>;</span><br><span class="line">        N = (<span class="built_in">length</span>(pulses)+<span class="number">1</span>)/<span class="number">2</span>;</span><br><span class="line"></span><br><span class="line">    <span class="keyword">elseif</span> pulses(<span class="number">1</span>) == <span class="number">-1</span> &amp;&amp; pulses(<span class="keyword">end</span>) == <span class="number">1</span></span><br><span class="line">        t2 = <span class="number">0</span>;</span><br><span class="line">        t3 = <span class="number">0</span>;</span><br><span class="line">        N = (<span class="built_in">length</span>(pulses)+<span class="number">2</span>)/<span class="number">2</span>;</span><br><span class="line"></span><br><span class="line">    <span class="keyword">elseif</span> pulses(<span class="number">1</span>) == <span class="number">1</span> &amp;&amp; pulses(<span class="keyword">end</span>) == <span class="number">-1</span></span><br><span class="line">        t2 = index(<span class="number">1</span>)/fs;</span><br><span class="line">        t3 = (wl-index(<span class="keyword">end</span>))/fs;</span><br><span class="line">        N = (<span class="built_in">length</span>(pulses))/<span class="number">2</span>;</span><br><span class="line"></span><br><span class="line">    <span class="keyword">else</span></span><br><span class="line">        <span class="built_in">disp</span>(<span class="string">&#x27;This should not be happening!&#x27;</span>)</span><br><span class="line">        keyboard; <span class="comment">%better to debug</span></span><br><span class="line">    <span class="keyword">end</span></span><br><span class="line">       </span><br><span class="line">    s4 = <span class="built_in">find</span>(aux(<span class="number">2</span>*wl+<span class="number">1</span>:<span class="number">3</span>*wl)==<span class="number">1</span>);</span><br><span class="line">    <span class="keyword">if</span> <span class="built_in">isempty</span>(s4)</span><br><span class="line">        t4 = <span class="number">1</span>;</span><br><span class="line">    <span class="keyword">else</span></span><br><span class="line">        t4 = s4(<span class="number">1</span>)/fs;</span><br><span class="line">    <span class="keyword">end</span></span><br><span class="line">        </span><br><span class="line">    tci6(<span class="built_in">j</span>+<span class="number">1</span>) = <span class="number">1000</span>/((N<span class="number">-1</span>)+(t2/(t1+t2))+(t3/(t3+t4)));</span><br><span class="line">    </span><br><span class="line">    <span class="comment">%Plot data</span></span><br><span class="line">    <span class="keyword">if</span> verbose</span><br><span class="line">        </span><br><span class="line">        f = <span class="built_in">figure</span>;</span><br><span class="line">        t = [wsamples1 wsamples2 wsamples3]/fs;</span><br><span class="line">        stage = [stage1&#x27; stage2&#x27; stage3&#x27;];</span><br><span class="line">        </span><br><span class="line">        subplot(<span class="number">211</span>)        </span><br><span class="line">        <span class="built_in">plot</span>(t,stage); <span class="built_in">hold</span> on;</span><br><span class="line">        <span class="built_in">plot</span>(t,[th1*<span class="built_in">ones</span>(<span class="number">1</span>,wl) th2*<span class="built_in">ones</span>(<span class="number">1</span>,wl) th3*<span class="built_in">ones</span>(<span class="number">1</span>,wl)],<span class="string">&#x27;r&#x27;</span>);</span><br><span class="line">        xlabel(<span class="string">&#x27;time/s&#x27;</span>);</span><br><span class="line">        ylabel(<span class="string">&#x27;ECG and threshold&#x27;</span>);</span><br><span class="line">        </span><br><span class="line">        subplot(<span class="number">212</span>);</span><br><span class="line">        <span class="built_in">plot</span>(t,becg,<span class="string">&#x27;k&#x27;</span>); <span class="built_in">hold</span> on; stem(t,aux,<span class="string">&#x27;r&#x27;</span>); <span class="built_in">hold</span> on;</span><br><span class="line">        line([<span class="built_in">j</span>+<span class="number">1</span> <span class="built_in">j</span>+<span class="number">1</span>],[<span class="number">-1.2</span> <span class="number">1.2</span>]);line([<span class="built_in">j</span>+<span class="number">2</span> <span class="built_in">j</span>+<span class="number">2</span>],[<span class="number">-1.2</span> <span class="number">1.2</span>]);</span><br><span class="line">        axis([t(<span class="number">1</span>) t(<span class="keyword">end</span>) <span class="number">-1.2</span> <span class="number">1.2</span>])</span><br><span class="line">        xlabel(<span class="string">&#x27;time/s&#x27;</span>);</span><br><span class="line">        ylabel(<span class="string">&#x27;pulses&#x27;</span>);</span><br><span class="line">        msg = sprintf(<span class="string">&#x27;t_1=%2.2f\t\t t_2=%2.2f\t\t t_3=%2.2f\t\t t_4=%2.2f&#x27;</span>,...</span><br><span class="line">            t1,t2,t3,t4);</span><br><span class="line">        title(msg)</span><br><span class="line">        text(<span class="built_in">j</span>+<span class="number">1.2</span>,<span class="number">-0.5</span>,[<span class="string">&#x27;TCI = &#x27;</span> num2str(tci6(<span class="built_in">j</span>+<span class="number">1</span>))])</span><br><span class="line">        </span><br><span class="line">        </span><br><span class="line">        <span class="built_in">hold</span> off</span><br><span class="line">        <span class="comment">% 是否进行调试，观看效果</span></span><br><span class="line">        keyboard;</span><br><span class="line">        close(f);</span><br><span class="line">    <span class="keyword">end</span></span><br><span class="line"></span><br><span class="line">    becg1 = <span class="built_in">zeros</span>(<span class="number">1</span>,wa/<span class="number">3</span>);</span><br><span class="line">    becg2 = <span class="built_in">zeros</span>(<span class="number">1</span>,wa/<span class="number">3</span>);</span><br><span class="line">    becg3 = <span class="built_in">zeros</span>(<span class="number">1</span>,wa/<span class="number">3</span>);</span><br><span class="line">    </span><br><span class="line"><span class="keyword">end</span></span><br><span class="line"></span><br><span class="line">tci = <span class="built_in">mean</span>(tci6);</span><br></pre></td></tr></table></figure>

<p>TCI算法</p>
<p><img src="H:\myboke\mybike\source\images\TCI算法得到的图" alt="image-20200520195855735"></p>
<p>准确率为100%结束</p>
<p><img src="http://img1.cache.netease.com/96/2014/8/12/2014081209202754754.jpg"></p>

    </div>
</article>


                </div>
                <aside class="col-md-4 gal-left" id="sidebar">
    <!-- 此为sidebar的搜索框, 非搜索结果页面 -->
<aside id="sidebar-search">
    <div class="search hidden-xs" data-aos="fade-up" data-aos-duration="2000">
        <form class="form-inline clearfix" id="search-form" method="get"
              action="/search/index.html">
            <input type="text" name="s" class="form-control" id="searchInput" placeholder="搜索文章~" autocomplete="off">
            <button class="btn btn-danger btn-gal" type="submit">
                <i class="fa fa-search"></i>
            </button>
        </form>
    </div>
</aside>
    <aside id="sidebar-author">
    <div class="panel panel-gal" data-aos="flip-right" data-aos-duration="3000">
        <div class="panel-heading" style="text-align: center">
            <i class="fa fa-quote-left"></i>
            esy
            <i class="fa fa-quote-right"></i>
        </div>
        <div class="author-panel text-center">
            <img src="/imgs/avatar.jpg" width="140" height="140"
                 alt="个人头像" class="author-image">
            <p class="author-description"><p>esy</p>
</p>
        </div>
    </div>
</aside>
    
    <aside id="sidebar-recent_comments">
    <div class="panel panel-gal recent hidden-xs" data-aos="fade-up" data-aos-duration="2000">
        <div class="panel-heading">
            <i class="fa fa-comments"></i>
            最新评论
            <i class="fa fa-times-circle panel-remove"></i>
            <i class="fa fa-chevron-circle-up panel-toggle"></i>
        </div>
        <ul class="list-group list-group-flush"></ul>
    </div>
</aside>
    
    <!-- 要配置好leancloud才能开启此小工具 -->
    
    
    <aside id="sidebar-recent_posts">
    <div class="panel panel-gal recent hidden-xs" data-aos="fade-up" data-aos-duration="2000">
        <div class="panel-heading">
            <i class="fa fa-refresh"></i>
            近期文章
            <i class="fa fa-times-circle panel-remove"></i>
            <i class="fa fa-chevron-circle-up panel-toggle"></i>
        </div>
        <ul class="list-group list-group-flush">
            
            <li class="list-group-item">
                <span class="post-title">
                    <a href="/2020/11/05/python%20work/%E6%9C%80%E5%B0%8F%E4%BA%8C%E4%B9%98%E6%B3%95/">最小二乘法</a>
                </span>
            </li>
            
            <li class="list-group-item">
                <span class="post-title">
                    <a href="/2020/11/05/python%20work/%E7%BB%9F%E8%AE%A1%E5%AD%A6%E4%B9%A0-%E7%AC%AC%E4%B8%80%E7%AB%A0/">统计学习--第一章</a>
                </span>
            </li>
            
            <li class="list-group-item">
                <span class="post-title">
                    <a href="/2020/11/04/python%20work/hello-world/">Hello World</a>
                </span>
            </li>
            
            <li class="list-group-item">
                <span class="post-title">
                    <a href="/2020/11/03/python%20work/%E5%88%86%E7%B1%BB%E6%A8%A1%E5%9E%8B%E7%9A%84%E8%AF%84%E4%BC%B0%E6%8C%87%E6%A0%87/">分类模型的评估指标</a>
                </span>
            </li>
            
            <li class="list-group-item">
                <span class="post-title">
                    <a href="/2020/10/21/python%20work/10-21-%E7%88%AC%E8%99%AB%E5%9F%BA%E7%A1%80/">10-21 爬虫基础</a>
                </span>
            </li>
            
            <li class="list-group-item">
                <span class="post-title">
                    <a href="/2020/09/25/python%20work/CRF%E7%9A%84%E6%95%B4%E4%BD%93%E6%B5%81%E7%A8%8B%E7%BB%93%E6%9E%9C/">CRF的整体流程结果</a>
                </span>
            </li>
            
            <li class="list-group-item">
                <span class="post-title">
                    <a href="/2020/09/25/python%20work/nlp-crf%E6%A8%A1%E5%9E%8B/">nlp_crf模型</a>
                </span>
            </li>
            
            <li class="list-group-item">
                <span class="post-title">
                    <a href="/2020/09/25/python%20work/%E6%95%B0%E5%AD%A6%E5%BB%BA%E6%A8%A1%E9%97%AE%E9%A2%983/">数学建模问题3</a>
                </span>
            </li>
            
            <li class="list-group-item">
                <span class="post-title">
                    <a href="/2020/09/25/python%20work/%E6%95%B0%E5%AD%A6%E5%BB%BA%E6%A8%A1%E9%97%AE%E9%A2%982/">数学建模问题2</a>
                </span>
            </li>
            
            <li class="list-group-item">
                <span class="post-title">
                    <a href="/2020/09/25/python%20work/%E6%95%B0%E5%AD%A6%E5%BB%BA%E6%A8%A1%E9%97%AE%E9%A2%981/">数学建模问题1</a>
                </span>
            </li>
            
        </ul>
    </div>
</aside>
    
    
    <aside id="sidebar-rand_posts">
    <div class="panel panel-gal recent hidden-xs" data-aos="fade-up" data-aos-duration="2000">
        <div class="panel-heading">
            <i class="fa fa-refresh"></i>
            随机文章
            <i class="fa fa-times-circle panel-remove"></i>
            <i class="fa fa-chevron-circle-up panel-toggle"></i>
        </div>
        <ul class="list-group list-group-flush">
            
            
            <li class="list-group-item">
                <span class="post-title">
                    <a href="/2020/07/04/python%20work/HMM%E6%A8%A1%E5%9E%8B%E7%9A%84%E5%AD%A6%E4%B9%A0/">HMM模型的学习</a>
                </span>
            </li>
            
            <li class="list-group-item">
                <span class="post-title">
                    <a href="/2019/12/10/python%20work/matplotlib%E5%9B%BE%E5%83%8F%E7%BB%98%E5%88%B6%E5%AD%A6%E4%B9%A0/">matplotlib图像绘制学习</a>
                </span>
            </li>
            
            <li class="list-group-item">
                <span class="post-title">
                    <a href="/2020/09/25/python%20work/nlp-crf%E6%A8%A1%E5%9E%8B/">nlp_crf模型</a>
                </span>
            </li>
            
            <li class="list-group-item">
                <span class="post-title">
                    <a href="/2020/02/15/python%20work/python-wfdb%E4%BD%BF%E7%94%A8/">python_wfdb使用</a>
                </span>
            </li>
            
            <li class="list-group-item">
                <span class="post-title">
                    <a href="/2020/03/25/python%20work/python%E2%80%94%E2%80%94%E8%AF%BB%E5%8F%96excel%E6%96%87%E4%BB%B6%E4%B9%8B%E5%A5%87%E8%91%A9%E6%95%B0%E6%8D%AE/">python——读取excel文件之奇葩数据</a>
                </span>
            </li>
            
            <li class="list-group-item">
                <span class="post-title">
                    <a href="/2020/02/09/python%20work/python%E6%A8%A1%E5%9D%97%E7%9A%84%E6%9E%84%E5%BB%BA%E5%92%8C%E8%B0%83%E7%94%A8/">python模块的构建和调用</a>
                </span>
            </li>
            
            <li class="list-group-item">
                <span class="post-title">
                    <a href="/2020/09/02/python%20work/tensorflow2-3%E7%89%88%E6%9C%AC%E5%AE%89%E8%A3%85/">tensorflow2.3版本安装</a>
                </span>
            </li>
            
            <li class="list-group-item">
                <span class="post-title">
                    <a href="/2020/08/08/python%20work/%E5%B0%86python%E8%BD%AC%E6%8D%A2%E4%B8%BAexe/">将python转换为exe</a>
                </span>
            </li>
            
            <li class="list-group-item">
                <span class="post-title">
                    <a href="/2019/12/16/python%20work/%E5%BF%83%E7%8E%87%E5%8F%98%E5%BC%82%E6%80%A7/">心率变异性</a>
                </span>
            </li>
            
            <li class="list-group-item">
                <span class="post-title">
                    <a href="/2020/09/07/python%20work/%E6%95%B0%E6%8D%AE%E9%9B%86%E4%BB%8B%E7%BB%8D/">数据集介绍</a>
                </span>
            </li>
            
        </ul>
    </div>
</aside>
    
    
    <aside id="gal-sets">
        <div class="panel panel-gal hidden-xs" data-aos="fade-up" data-aos-duration="2000">
            <ul class="nav nav-pills pills-gal">

                
                <li>
                    <a href="/2020/05/18/python%20work/VF%E5%9F%BA%E7%A1%80%E7%9F%A5%E8%AF%86/index.html#sidebar-tags" data-toggle="tab" id="tags-tab">热门标签</a>
                </li>
                
                
                <li>
                    <a href="/2020/05/18/python%20work/VF%E5%9F%BA%E7%A1%80%E7%9F%A5%E8%AF%86/index.html#sidebar-friend-links" data-toggle="tab" id="friend-links-tab">友情链接</a>
                </li>
                
                
                <li>
                    <a href="/2020/05/18/python%20work/VF%E5%9F%BA%E7%A1%80%E7%9F%A5%E8%AF%86/index.html#sidebar-links" data-toggle="tab" id="links-tab">个人链接</a>
                </li>
                
            </ul>
            <div class="tab-content">
                
                <div class="cloud-tags tab-pane nav bs-sidenav fade" id="sidebar-tags">
    
    <a href="/tags/py-study/" style="font-size: 18.446350830709px;" class="tag-cloud-link">py_study</a>
    
    <a href="/tags/nlp/" style="font-size: 19.72672041823357px;" class="tag-cloud-link">nlp</a>
    
    <a href="/tags/Graduation-work/" style="font-size: 13.519227005397843px;" class="tag-cloud-link">Graduation work</a>
    
    <a href="/tags/work/" style="font-size: 19.021384801021906px;" class="tag-cloud-link">work</a>
    
    <a href="/tags/hexo/" style="font-size: 11.139277258985548px;" class="tag-cloud-link">hexo</a>
    
    <a href="/tags/%E4%B8%AA%E4%BA%BA%E5%8D%9A%E5%AE%A2%E6%90%AD%E5%BB%BA/" style="font-size: 10.783683165088123px;" class="tag-cloud-link">-个人博客搭建</a>
    
    <a href="/tags/malab-%E6%AF%95%E4%B8%9A/" style="font-size: 13.230057362665415px;" class="tag-cloud-link">-malab -毕业</a>
    
    <a href="/tags/python-%E4%BA%BA%E5%B7%A5%E6%99%BA%E8%83%BD/" style="font-size: 15.7041296102277px;" class="tag-cloud-link">-python -人工智能</a>
    
    <a href="/tags/python/" style="font-size: 16.693593666406645px;" class="tag-cloud-link">python</a>
    
    <a href="/tags/python/" style="font-size: 16.276129262353226px;" class="tag-cloud-link">-python</a>
    
    <a href="/tags/mathematical-modeling/" style="font-size: 11.541334523563243px;" class="tag-cloud-link">mathematical modeling</a>
    
    <a href="/tags/statistical-learning/" style="font-size: 11.753624441482522px;" class="tag-cloud-link">statistical learning</a>
    
</div>
                
                
                <div class="friend-links tab-pane nav bs-sidenav fade" id="sidebar-friend-links">
    
    <li>
        <a href="http://kdays.net/days/" target="_blank">KDays Forum</a>
    </li>
    
    <li>
        <a href="http://www.gal123.com/" target="_blank">绅士导航♂</a>
    </li>
    
    <li>
        <a href="http://www.moe123.com/" target="_blank">萌导航</a>
    </li>
    
</div>
                
                
                <div class="links tab-pane nav bs-sidenav fade" id="sidebar-links">
    
    <li>
        <a href="https://github.com/ZEROKISEKI/" target="_blank">Github</a>
    </li>
    
    <li>
        <a href="https://coding.net/u/SORA1" target="_blank">Coding</a>
    </li>
    
    <li>
        <a href="https://www.zhihu.com/people/aonosora/activities" target="_blank">知乎</a>
    </li>
    
</div>
                
            </div>
        </div>
    </aside>
    
</aside>
            </div>
        </div>
    </div>
    <footer id="gal-footer">
    <div class="container">
        Copyright © 2018 esy Powered by <a href="https://hexo.io/" target="_blank">Hexo</a>.&nbsp;Theme by <a href="https://github.com/ZEROKISEKI" target="_blank">AONOSORA</a>
    </div>
</footer>

<!-- 回到顶端 -->
<div id="gal-gotop">
    <i class="fa fa-angle-up"></i>
</div>
</body>

<script src="/js/activate-power-mode.js"></script>

<script>

    // 配置highslide
	hs.graphicsDir = '/js/highslide/graphics/'
    hs.outlineType = "rounded-white";
    hs.dimmingOpacity = 0.8;
    hs.outlineWhileAnimating = true;
    hs.showCredits = false;
    hs.captionEval = "this.thumb.alt";
    hs.numberPosition = "caption";
    hs.align = "center";
    hs.transitions = ["expand", "crossfade"];
    hs.lang.number = '共%2张图, 当前是第%1张';
    hs.addSlideshow({
      interval: 5000,
      repeat: true,
      useControls: true,
      fixedControls: "fit",
      overlayOptions: {
        opacity: 0.75,
        position: "bottom center",
        hideOnMouseOut: true
      }
    })

    // 初始化aos
    AOS.init({
      duration: 1000,
      delay: 0,
      easing: 'ease-out-back'
    });

</script>
<script>
	POWERMODE.colorful = 'true';    // make power mode colorful
	POWERMODE.shake = 'true';       // turn off shake
	// TODO 这里根据具体情况修改
	document.body.addEventListener('input', POWERMODE);
</script>
<script>
    window.slideConfig = {
      prefix: '/imgs/slide/background',
      ext: 'jpg',
      maxCount: '6'
    }
</script>

<script src="/js/hs.js"></script>
<script src="/js/blog.js"></script>



<script src="/js/oni.js"></script>




</html>